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ABSTRACT 


A simple accurate weight- function W(m.) has been 
used instead of the exact wOx) and with this weight- 
function orthogonal polynomials have been developed. 
The solution of the one-speed transport equation has 
been found out by talcing the transient integral as a sum 
and using a rational function approximation (P^(VfU) 
to the exact This method of finding out the 

solution of the transport equation by approximation 
has been applied to different problems of the transport 
theory and the results are found to be reliable. 
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CHAPTER 1 


IHTRODUCTIDN 


The one-speed transport equation plays an important 
role in the transport theory. The complete set of eigen- 
functions of this equation contains singular distributions 

in the form of the Cauchy principal value and dirac delta 
<jki sfevibv-ttorii 

fun& fel^ nals besides regular functions, see Case and 
Zweifel £l] for details. 


We are concerned with the solution of the hcmogeneous 
one-speed equation 


3x 


+ ?(x,u.) = I 


1 

-1 


(i.i) 


for full-and half-range boundary conditions. There are 
two different ways of obtaining the solution of this equation. 
One is the Case's method and the other is Wiener-Hcpf 
method. The Case's method, like the wiener-Hbpf method, 

eaCOiCfc 

is aeeurat-e and gives the exact solution of Eq, Cl.l;. 

It is analogous in some respects to the method of separation 
of variables ccamonly used for the solution of partial 
differential equations. In both, a complete set of elementary 
basis solutions is sou^t, and then a suitable cxxnbination 
of solutions is found that will satisfy the boundary 
conditions or the conditions at the source. The major 
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difference is that one of the elementary solutions of the 
transport equation is singular, i*e., it has meaning only 
when it appears inside integrals- 

The Case's solution in standard notation [l] is 



®ie necessary and sufficient conditions for the determi- 
nation of the constants a^^ and aCv) are the orthogonality 
integrals [l] 



3 


/ jj. (p^ (ul}611 = N 

o± ot 


(jLL)cJl = N(v)6(.^V') V,V*€(^1,1> 


where 


cv 


N 


0 + 


i 2 


o C • 


2^0-1 

o 




2> 


N(a>) = y[(l- |2i In ] 


for the full-range problem, and when VfV' e (o,l). 


/ w ( m .) (P ,(fi)dli = O 

o o+ V 


/ W^) = w(y) (N(v)/a))5Cy-a)') 


1 2 
/ w(fi.) <P^^O^)dJ, = + Cf cv^) x(+ Vq) 

O "• 


for tlffi half~raftge case fl]« 

The set of eigen-functions —1 ■* 

of Eq»(l«l) are ccanplete in L^j^OJl] with respect to 
weight function wOJ-) given by 

= sftc) ■nr^TTiTOT 

where 

x(-p) SB 


tl.6a) 

(i.6h) 

(1.6c) 

(l.6d) 

(1.7a) 

(1.7b) 

(1.7c) 

y <11 

the 

( 1 - 8 ) 

(1.9) 
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cv 1 -1 

Q (-11 ) =1- ju. / — 

^ o Oi+t)a(*.t) 


dt (1.10) 


with 


X(0) = ^ 




( 1 . 11 ) 


Also 


vfCl-c)-! 

Xiv) =-( ^ 


2a^a-c) 


1/2 
■) ' 


( 1 -. 12 ) 


x{-v ) = xCo)^) 

O OM O 


(1-13) 


where 


®oM = -e=®(-2V’'o> ' 


(1.14) 


being the extrapolated and point- Once the coefficients 
are determined^ then the flux and leakage can be calculated 
according to the problem considered. The weight function 
W(]i) given by Eq, (1.8) is a non-analytical function* 

Also tl^ weight function is expressed in terms of fi(-4i.) 
which satisfies the non-linear integral equation (l-lQ), 

Though the Case's method solution is exact# it is quite 
complicated to evaluate the integral in the solution of the 
transport equation because of its singular nature and the 
nature of aCj?). For exanple# consider the Milne problaa £l] 
in which we seek solutions of the bOTOgenpus transport 
equation stibjeGt to 
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■* # X** «® (l»i5) 

vdiere f(x,/J-) Is the solution of the problem, and to 
the condition at x = O, 

f (0,JJ.) =0 , )J. > 0 ■ (1.16) 


The solution to this problem can be taJcen to be a linear 
combination of the fundanental solutions which vanish at 
infinity plue . Hence,' 

1 

I^Cx/W-) = ^ (x,M-) '+ a 4.^-. (x,»i> + / 'A(a>)^,,(x,fi)d» 

O** Ot Oj* ^ y 

(1.17) 

The condition (1.16) then gives 

1 

(m>) = a„, + / A(5;)<P^,(fi)d3; , ix>0 (1.18) 

O— 0+ 0+ V. ^ mm 

O 

By applying the orttogonality relations [l], and Aiv) 

-can be jBound Immediatelyi 

a(i/) = - / w(M-) <P^_(|i.)dM. (1.19a) 

o 

or 


2 

X(-y ) <p (j>)3; 

A(y) = —2 

(y_-y)y{y)N(y) 

o 

where 

2 

nCv) s y[ (l-.ca» tan h“^y> + c^n^v^/4] 


Cl.l9b) 


(1.19c) 
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Similarly 

,2 1 . 

“[ f = - / Y <*>^,0^)0^ 

(1.20a) 

so that 

= x(-v )/X(v ) {l.20b) 

In the Albedo problem fl]^ in which we wish to obtain 
a solution to the homogeneous transport equation dBor 
O <x subject to the twD boundary conditions 

f^{0,iu) = 6 0^-^^) > O Al.21a) 

and 

lim W (x,h) - O, (1.21b) 

X “ 

the solution is taken in the form 

1 

5'^(x,j^) a + / A(v)W^(x^|Lt)d3; , (1.22) 

o 


vhere, the coefficients a^_j_ and a(v) are obtained from 
the boundary conditions, Eq. (l.21a) using the half -range 
orthogonality formulas [l]. 


VWQJ^) 

ACy) a 2L-0_ 

N(y)w(a?) 

and 

= -2y {M._)/cy^ x(v ) 


(1»23> 

(1424) 



Thus from the two examples mentioned ahove, ws see 
that the nature of Aiv) aee quite ccmpllcated and the 
iiitegrateen involving this In Eqs. (l«17) and Cl»22) 
camot be evaluated analytically. Therefore It beccines 
necessary to find an approximation for the integral in 
the solution of the one-speed transport equation. This 
Is the principal objective of the present thesis. 

In the conventional approximation, the exact 

solution ^(x,h) given by Eqi- (1.2a) is approximated in 
the form of a truncated Legendre polynomial expansion of 
the functions and and this forms the 

O-r V ^ 

basis of the approximation [ 2 ]. This approximation 

is equivalent to 


gjj+iCv) = 0, (1.25) 

where Ihe coefficients in the expansion of 

?(X/h) as a Sim and are given by 

1 

g (y) = r P„(h)dlLi (1.25a> 

il y H 

Using Eq. (1.3 ) and (1.4 ) for we have 



(1.27) 



itihere Legendre polyixxnials of secox^ Icind [ 2 ] • 

The Eq, (1.25) discretizes it to (N+l) values* The 
conventional approximation then considers the largest 

tv)0 of them as the asymptotic relaxation lengths and the 
remaining (N-l) as the transient ones. The two largest 
roots are obtained from Eq, (1,26) by imposing the condition 
that g^Cv) =1* 

In the transport- theoretic P.-(TF ) approximation, 

JM jj . 

the infinite sums for the two asymptotic roots, i,e*, 

. (ji) are retained and the series for the transient part 
04 ; 

is terminated. The solution is varitten as [^2] 

= a e 3 {p(-x/i> ) (W’) + a_ expix/v) 9 (m-) 

^ 0+ ' o 0+ o— ^ o o— 

N 

+ 2 aCv .)exp(-3(/yj) 2 g (y } P On ) 

J=roots ^ J n=0 2 ''n J n 

( 1 , 2 &) 

The coefficients in Eq, (1*28) are to be determined from 
the boundary conditions. Here the standard free surface 
Marshalc boutKiary condition is modified to get the exact 
extrapolated end point* The solution of the soxirce free 
Milne problem [2} is considered and the integral in the 
solution is replaced by a sum and the coefficients are 
obtained from the half-rai^e orthogonality relations {^2]. 

This TP^ formalism has be ^ applied bo the following 
problems £ 2 ] 
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The emergent, angular distribution and le^age for 
the Milne problem 
2s The scalar flux in Milne problem when c = 1 

3» leakage in cx>nstant source Milne problem 

4« the critical problem 

5. Interface and source problems 

and the results are ccmpared with the exact values* The - 
TPjj approximation gives significant improvements over 
conventional method# except in those cases where the 
effect of the transients is especially importeait. 

Anoth^ method known as the Fj^ method in Neutron- 
Transport theory [6] is used to solve the half-space 
albedo problem and the half-space constant-source Milne 
problon. The half-space problem is defined# for c < 1# 
by 

ja + ¥Cx#fx) = f / i'(x#uOoH' + a Cl.29> 

W(0,M.> := 1-a ii > O (1*30) 

and 

^(x#M ) •* # as X (l*3l) 

Here a sa 1 yields the usual constant-^urce Milne problem# 
and assO yields the half-space albedo i^oblem» A particular 
solution of the problem is £73 



1 

W<x,fi) as A{v^)<P(5;^,M)e3<p(-3c/v^) + / A(y)<l>(v,4)«3^(-3^y)d» 

o 

+ ^ , H > O (1*32) 

Then, on noting that 
1 

/ [f {3C,;J.) - JA <P (-r,|i)ai =0, s e p (1».33> 

-1 

where i e P «> t or 5 « v e (o, l), we have 

^ 9 (5,M-)^(o,4i)P' <ii. = + {l-a)£l-S log(l+ ^>3 

I e P Cl.34) 

The last equation has been obtained splitting the interval 
(-1,1) to (-1,0) and (O, l) and using the Eq. (1.30) in the 
last integral, ij.e., for the interval (0,l), 

We enter for the emergent distribution, 

N 

5'C0,-m) = IT , )UL > O (1.35) 

<X=D 

into Eq, ( 1 , 34 ) evaluated at selected values of S to obtain 
the Pjj equations 

^ 9= V 

+ (W)[l-tp log<l+ ^1] / 

^= 0,1,2/. ../N (1,36) 


%^ere 



(I*36a) 


with 

= ~ -1-1 log(l+ ^) St log -1 (l*36b) 

o 

Then * v^, = o, = i# and the remaining tp are 

spaced equally in the interval fO,l]j and is or» possible 
set of collocation points that has been used in the 
literature [6]. 

The leakage given by 

1 

j B / fxwco, 

-1 

6 1 
s / jULW(0,ji.)cM. + / txT(0,fi)^ 

-1 o 

1 

* / )xfCo,-ju)dA (l,37> 

o 

is calculated for the cases 0=0 and a « 1 and oora^red 
with the exact values* This method/ although paiHdcularly 
concise, yields excellent ntamerical results for the prol:d.ems 
oonsid^ed* 

In the present thesis, we utilise tl» coiKJept of a 
rational function approximation introduced by Sengupta [ 3 ], 
[43 for the singular eigen-distribution fo reduce the 
transient integrals of Eq. {l*2b) to a sum. The Itegendre 
or Chebyshev polynomials are usiis^y used as the complete 
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set and their zeros as collocation points for obtaining the 
approximate solutions of Neutron Transport equation* Since 
the actual set consists of the case eigen— functions w^^r^t* 
the weight function wOj.) it is rx>t the best possible to 
use the classical, orthogonatl polynomials in transport theory*. 
Therefore it becones necessary to orthogonal is e the 
fundamental set w.r.t* wCut) and use this as a 

proper substitute for the case eigen-f unctions and 

[ 5 ]. 

In our present work, a weight function is used 

in the calculations instead of given by Ekj* (1*8) to 

obtcdn the coefficients occuring in the solution of the 
transport equation and this has been properly justlf i^ in 
the next chapter. A new set of orthogonal polynaaials 

w*rvt* the weight function has been constructed 

Ref . 1^5 J and the values of M and v are tak^ as the zeros 
of these polynomials depending on the niaaber of traz^ient 
terms considered* Then the coeffici®:its are obtained for 
the problems given below 

(i) standard Milne problem 

(ii) constant-source Milne problem 

(iii) half-space Albedo problem 

mad then the emergent angular distribution and le^e^e are 
calailated and ccmpared with the esact values. This method 



chapter: 2 


METHOD formulation 

1* TRANSPORT THEORETIC APPRO XIMATIOH [3], [ 4 ], [s], [S] 

The solution of the one-speed equation |33 is given by 
Eqs« {l*2a) and (l*.2b) and the Case eigen- functions are 
given by Eqs- (l«3) and The constants a^^^ and 

a(i>) occuring in the solution are determined by using the 
orthogonality integrals given by Eqs. (l«6) and (l.7> 
for the full-and half-range problens. Here X-furK^tion is 
given by Eq» (l«9) and W- function by Eq, (1*8), 

In this approximation, it is reasoned as follows [3] 

(i) The discretization of V e (-1,1) or (0,l) must be 

done such that these roots, a>j,j=l, 2... are consistent 
with the parti ciilar choice of The P^^ approximation 

has this desirable property, as all its (N+l) Vj are 
soluti QTiR of a polynomial equation of degree N+l, Eq* Cl»2S) 
In contrast, by requiring the asynptotic to satisfy 

Eq. Cl«5) and the transient the equation (1*25), the 
TPj^ procedure £2^ violates this property* 

(ii) The natural basis functions for the solution of the 
transport Equation are the case eigen- functions 

(fi)," such as the fPjj(P'))j^^ 





(as in the case)> or a ccHtibination 

CPjjOj-)}} (as in the TPj^ case), is likely to be unsatisfactory* 

(iii) By the very nature of the eigen- functions, a rational 
function approximation is expected to be superior to a 
standard polyncmial approximation- 

In the case of a half -range problem, a single but 
remarkably accurate expression for w0j.),[3] 


wOi.) 


c 

“ 2Q ( 1-c) 


Uv V l-G + 

^ O ’ 


( 2 . 1 ) 


where 0 is taken to be a constant has been constructed [3] 
for use in the orthogonality integrals (1*7), The value 
of Q is obtained by evaluating the _ orthogonality integral 
(l#7c) • Also with this constant Q, a first iterant of 
Q(-U) from Eg. (l.lO)-and using this, x(-u) from Eg,. (1.9)- 
are obtained [3]. 


Now, the solution of the transport e{guation (l.l) can 
be written as 


^(x,U) = exp(^x/v^) 

+ S A(y.) exp(— 39^^.) ^g(a’4,u) (2#2) 

where ^p (y,Ji-) is the proper rational function approximation 
to ^P(v,M’), V £ (-1,1) and is given by [4] 



£ 


IP (vm) = — indil L + 2is<y.>. , 

^ ' 2 (1.^)2+ ^2 "e [(WP)2«2j' 

V e (-1,1) (2,3) 

The constants in Eg. (2.2) can be determined using the 
boundary conditions and orthogonality relations satisfied 
by the set { The function <P^(y,|Ji-), 

as e •• 0 , tends to the sum of the tw distributions that 
qomprise One is the ftrincipal value distribution 

and the other is the delta distribution. 


Using the normalisation of <Pg/ we get TveCv) from 


Eg# ( 2 * 3 ), 

?v£(i>) = 


H, 


tan"^ + tan“^ 

e e 


[ 1 - 


Cl+v)^+e^ 


CV ■, Vi+y/ -t-t t 

4~ In — 72 — oJ 
^ (l-v)^+e^ 


( 2 . 4 ) 


where 

Tte = 2 tan”^ 1/e (2%5) 

and Ae(i>) •* X(y) as e -• 0 . 

The convergence of the tv30 distributions that ccanprlse 
(a>,]Li) has been verified [S] by finding the solution of 
love's Integral eguation 

1 U^(t) 

U^(s) = 1 - - / — 5 - dt (2.6) 

^ 1“ -1 (s-t)2+ e2 

and the Cauchy Principal integral 



dt = 

o 


/ 


(2.73 



Mtl 

t-x 


where is a constant, 

o 

In Qiir present work, tl^ solution of the one-speed 
equation (l.l) is taken in the form given by Eq. (2.2). 

Then the given boundary condition is made use of and the 
coefficients and {aCvj)} ; are determined by 

considering the fJ-j's as the roots of Cjj_^^(ia^) r= O and 
Vj's as the aroots of Cjj(a;j) = O where Cj^(M-) are the 
orthogonal polynomials constructed in [s] and details of 
v^Oh are given in next section. The value of e is 
determined from the relation given by Eq, 

e = (2.7a) 

2N tan“-^ 1/e 

depending on the number of transient terms used* Once 
the coefficients are determined then the emergent angular 
distribution and Leakage can be calculated for various 
problems considered. 


2. JUSTIFICATION OF WEIGHT FUNCTION 

A weight function w(ii) given by Eq, (2.1) has been 
cxjnstructed in [3]. Considering the simplicity of w(m-) as 
cx»npared to the exact wOi) given by Eq, (1,8), the acchiracy 
of the former is remarkable'. With this w(^), the 
orthogonality integral (l.'i7c) are evaluated to give 



afflllc) % (<xc^+C2^ = -^1 


i2»8) 


and 


‘2Q (i-c) % (aD^+D^) = 


(fcy^) x(-a>^) 


(2.9) 


where. 


c, = 


C„ = 


D, = 


D„ = 


1 


2v (y -1) 
o o 

2cv ^ 
o 


I- 1 

- ^ 

2(y -l) 
o 

2cv^ 

o 

_1 

1 ■ r. 

2c3;“ 2v 

o o 


1 

-2— 4- 1 


V- 




2(y^+l) 2 c3^q 




The two equations (2.8) airf. (2:^9) are solved to get 
and 0^. I£ 0 is very nearly constant, which is assumed 
to be the case in Eq, (2.1), then the solution of the 
■above two equations for Q, i.e. and Q_, will be marly 

the same. To taJce care of the weak dependence of 0 on p. 
to a satisfactory degree, we use the average of 0 and 0 

Hr ■''••' 

i.e. let 

Q = i (£i_^ + 0 ) (2.10) 

As an, independent check on the utility of this w(p) in 
evaluating integrals of the type (1.7), the integral below 
was evaluated [3] as shown in Eq^jf (2.1l), 



1 , -1 
/ w(p.) [4U(l-c)] (OCF^^P^) C2*ll) 


v^ere 


P = / 
^ o 




n 


2 2 


dll 


1 

s= — C" T + + _ tr +•*•)+ i 

n— X n— 3 n— 5 


o 


n evBU 


v^ln[v^/ n odd 


and compared it with the exact value ^ ratio of the 

approximate to exact integrals for different c are 0.999594, 
0.999530, 0.999693 and 0.999832 for c = 0.^2, 0.4, 0.6, 0.8 
and 0.9 respectively. 

As a further check, we have evaluated the integral 1 
given by 


I = / w(ia) (P(v,li)di 

o 


(2.11a) 


and compared it with its exact value (for the real wQi)} 
of zero 


3 

c w 


I = 


80 (l-c)' ^1 4i) (l-c) (^^2 ^ ^2j 


[yv Vl-c + y^l 

o ^ J. (2.11b) 


where 


T dL In 4. B' -i- In + r.' In - 2 L 

“ 2 ^ ,,2 , 2v ^ v-1 + i^p 

“** jL V./ V 


(2.11c) 



A'' * - 


vjv '^1-c + vl 


B* s: - 




Cv?- v^) 


y[ y + V '/ 1--C ] 


(a;2_ j,2) 


The value of I has been obtained for different values of 
C and y and the outputs are attached in APPENDIX-A,- 
Considering the accuracy of these two examples, we can 
say tliat the simple analytic approximation C2.I) to the 
tr ans cendental weight function W (p ) , Eqs « (l. 8 )-Cl.lO), 
can be used with confidence. 

A first iterant of x(-M-) can be obtained from 
Eq* (1*9) with 0 (-M) replaced by the constant Q, i.e.> 


5^0 =551“ 


(2.12) 


A first iterant of OC-p) is obtained frcm Eq» (l.lO) as 


={ 1 -^ I W] 


(2*13) 


where 


^o*"^ I4i 

In -jj— 


W =* In 


+ oij in 



c 


a, = 


a. 


2y^Cyo-4i.) Cl-c) 


2a> iv -iiJ-) (l-c) 
o o 


a 


= x^(o) - 


(l-c)(y2^u2) 


witJi this a second iterant of xC-+l) is obtained as 




(2.14) 


The values of X^(-4i) given by Eq. (2.12)/ 3^(-jU.) 
given by Eq» (2.15) and the ERROR given by 


ERROR 


(EXACT VAIiUE - CALCiniiAT^ VALUE) x 100 
EXACT VALUE 


(2.15) 


are calculated for various values of ix and tl^ outputs are 
attached in APPENDIX We observe from the results that 
the error are very less and the above, therefore^ constitutes 
an extremely reliable approximation to Case's W and x 
functions . 

hm ^ 

3. GRTH3G0NAL POLYNCMIALS fs] 

£1/)X*4^} is orthogonalised w.r.t* W^) given by 
Eq* ( 2 . 4 )/ <1# using the prescription given by Golub 

and Weis eh* [ 9 3^. Let 

' 1^' ■' ^ . ■ . ' 

/ w(|^ ) p' ^ 3=1# 2# ♦ ♦ H-fl 


C 2 » 16 ) 



and 


where 


±-l 


’’ij = ^Ij - ’'vf i < i 


k=l 


ki 'kj' 


i-1 


’'11 = <“11 - \P 


2,V2 


Now construct 

i+1 


a. s 
i y 


'±^1,1 


i, i i-1 


( 2 . 17 ) 


( 2 . 18 ) 


^ » 

6 ~ 4tlf 1+1 ■}_•} O M 


i/i 


where y _= 1/ y_ , = O* 
0/0 O/i 


N 


(2.19) 


Then with C^=l/ C-1 = 0/ the set 


obtained from 


the recurrence relation 


Pi C^(M-) = (M—a^) ^=^'2/ .../N 

(2v20) 

form a complete orthonormal system in (0,1) w.r.t., w(M-)* 

For our present problem. 


M - „■■, c / ML 

" 2Q(l-cy [. y_-i^ ^ ^ 


(2>2l) 


o o 

which is exparessed in terms of 

1 

o o 






( 2 . 22 ) 


f'tl » ” ■ , ^ 


^ w • a««. I 


where 



O 


(2-23) 


^i+j-1 **■ ^i-j^ 


(2.24) 


The value of M . , can be evaluated recursively, ajod 

“‘"-J 

hence finally {C^Cm.)}^ is obtained frotn 

Eq. (2«20)fc The orthogonal polynomials C^(jJ.) idD,l, ..,,5" 
has been developed and presented below- 
By Eqs. (2.18), (2.19), (2.20) and (2.21), we have 
r, r. 


a - 1Ll 2 . p ^ 

“1 “ r, ' ^1 - 


• 1/1 


' 2,2 

'l/l 


r _ (M )V2 _ .1/2 
^ 1,1 " ^ 11 ^ " 


where 


5n = ^ ^ntl ® ^n..2 

^ i-C ^ 

^ = 2a'(X^^^ ' ^ 2a (i-c) 


and are given by Eq'.- (2,;^23) and (2v24) 


M 


12 


1,2 r, , ” ^ 1/2 

-L/J- Oq 


^ 2,2 = ^22 - ''1,2^ 

6? 1/2 

= ^«2 - 

o 


1/2 



“l “ ^ = 1“ ^^0^2“ 


1/2 


G^{m.) = Ox-.a^)/P^ 


Also 


= 


‘2,3 1,2 


2 r 


2,2 n,i 


« r- 

2 r 


‘ 3,3 

2,2 


^2,3“ ^23*'^1,2 ^1 3 ^'"^^2 2 


= (6_-r- « r, _)/r„ _ 

3 2 ■ X/ 3 ^ 2 jr 2 


^3,3* ^^33“^, 3 " ^2,3^ 


1/2 


7 ^ V2 

« (6^ — ^1/3 "* ^2# 3^ 


[Ox-a^) Cj^Cn) - ^jl/^2 
=: [ (iLL-a^) (ix-a^) - 


Similarly, ve have 


^ "^3,4 _ fXa 
3 ^ 3,3 ^ 2,2 


P. = 


^4,4 


3 T. 


3,3 




CgOi)* [(M-a^) C^O^} - §2 • 


“4 

1 
u 1 

[ 

_ ^3 ^ 4 

L"" 

1 

^3,3 


_ ^5,5 
^4,4 



^1,5 “ ^ 15 /^ 1 ,! " 

^2f 5 ~ ^25 "" ^1#2 ^ly5^/^2^2 
“ ^6g "* ^2.f2 ^2.f5^ ^^2^2 

^3^5 “ ^35 ” ^1,3 ^1,5 “ ^2,3 ^2#S^/^3,3 
"* ^^6 “ ^ 1,3 ^ 1,5 “ T 2^3 ^ 2 ^ 5 ^/^ 3,3 




Finally 


“5 

_ ^5 i6 __ 

” ^5^5 

ficS. . p _ 

^4,4 ^ 5,5 

^1>6 


,1 = 55 /So^'^^ 

^2,6 

= ‘«26 - 

^1,2 ^l,6^/^2,2 

^26 

= 5g 


^3,6 

= ^36 - 

r, „ r, , — r« « c}/'^r, 0 

1,3 1/0 2,3 2/0 ' 3/3 


= 67 



^ 4,6 “ ^ 46 *“^l /4 ^ 1,6 ' ^ 2,4 "^ 2 , 6 “ ^ 3,4 ^ 3 , 6 ^/^ 4,4 

%6 * % 





^1,5 

5-/6“ ^2,5 ^2,6“ 

^3^ 5^3, 6"^4, 5^4,^ 

^56 

69 




^ 6,6 ^ 

^66 - 

- 

2 2 
^2,6 •* ^3/6 ” • 

2 2 ^V2^ 

'^4,6“ ^5|.6 ^ 

M — » 

^66;. “ 

•^lo 




C^0i) = 


) C5OX) 

- P4 

» 


Note that 

C Oi) = 1, 
o 

The roots fof C^CM-) = 0 for N=sl to 5 have also been 
obtained by bisection method a»3 the program and the outputs 
are attached in APEENDIX-C. 

4. APPL ICATlOf^ 

The me-chod discussecl in section 1 of .this chapter 
has been applied to the following problems 

(i) Starxiard source- free Milne Problem £l] 

(ii) Con''tant-source Milne Problem [l] , 

(iii) . Half-space Albedo Problem [l] 
and they are discussed below- 

(i) Standar d source-fr ee Mi l m Problem 

The solution of the standard source-free Milne Problem 
is given by 



a expU30i.^> 

N 

+ S A(y.) exp(- 5 ^i> ) <P_(v.,u) (2.25) 

jsri -’ J t- J 

with boundary condition at x = O, 

W(o,u) = 0 , U> O (2*26) 

Note that in Eq* (2*25), the value of a^ is taken to be 
unity and 9^(v,IJ-) is given by Eq. (2.3). 

Then applying Eq, (2.26) to Eq, (2.25), we get 


%+‘*’o+^^ + A(a;.) <P^(.v.,ia) = o. (2,27) 

The coefficients a^, and Ihiv.)}^ , can be obtained hy 

o+ J J=1 

taking 4's as the roots of '-n+I^ = O, v*s as the roots 
of ^ the value of e frcxn Eq, (2#7a) for N 

equal to the number of transient terms considered. 


The emergent angular distribution If (0,4) is given 


W(0,4) = a^_j_ <Pj^_j,(4) +2 A(j)j)^^j (u), 

4> 0, y>0 (2,28) 

Note that for 4< O, W ** a® there is no 

chance of tto denominator of 



becoming zero. Hence the case of Principal value a»a Delta 
distribution doesn't arise* 



The leakage J is given by 
o 

J S5 / JJlf (0,JLt)(3J, 

-1 



o o 

/ u.<p (^.)dis + / I1<P (tJi.)dui + S aCv.) 
-1 -1 j J 


o 

/ I1<P (y.,^)m 
-1 ^ J 


cy 

04* 2 


J'^+l CV V 


cy. y .+1 

- S A(yj) [l-y^ In -^] 


(2.29) 


Note that in the last integral of the expression J, 
^g.Cyj#M-} - as ja < 0. 


(ii) Constant-source Milne -problem 

a) The solution of Eq, (l*29) with a5=l and boundary 
condition given by Eq, (l,30) and (1.31) is 




-a <p Ojl) exp(-x/y^) +S ACy. )9 (a;.//i)e>!pC-x/yj) 
o+ o+ ^ ' o . J E j j 



(2.30) 


Using the boundary conditions, the coefficients are obtained 
frcan the equation 

a <P (n) + S ACy J = ° (2.31) 

04- j ^ ^ ^ ^ 

as in the standard source— fir ee Milne j^roblem. 


After 
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determining the coefficients and { ACy^)} , the leakage 

is calculated by 

o 

J St / fi W(o,iJ.)dlJ. 

-1 



(2.32 > 

b) The solution of Eq?.. (1, 29) with a=0 and 
boundary condition given by Eqi, (l,30) and (l.3l) is 

expC-Vi'o^ ^ A(y.)<P^(y ,M.)e3^(-.3^yJ 

j J J J 

(2.33) 

using the boundary condition, the coefficients are 
detemined from 

a^^«P^^(h) +S A(y,) <3>_(j?.,p.) =®( (2.34) 

o+ o+ j J £ J 

as before* Then the leakage is obtained fran 
o 

J = / nT(0,IJi)<Si 
-1 

o o 

ss a J‘ fA'P (jJL)d|ti + S A(y,) J 

0+ 1 j ~1 ^3 
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5 . RESULTS AND CONCLUSION ' 

The values of and for various values of C 

ranging frcsn O.l(o.l) 0.9 are presented in Table 1» The 
value of Q is calculated from the Eq. ( 2 .I 0 ). 

The emergent angul^ distribution and the leakage 
for the standard source— free Milne problem are caloilated 
using Eq* (2*28) and (2.29) and the percentage errors 
using E<qv (2.15) and are presented In Tables 2 and 3» 

The exact values are t^en fran [ 2 ]. 

The leakage for the constant source Milne problem 
with a=D is calculated using Eq* (2*35) and the results 
with the percentage errors are presented in Table 4. The 

iw' 

exact values are taken frcm [7]. 

The leakage for the constant-source Milne problem 
with a=l is calculated using Eq, (2.32) and the results with 
the percentage errors are presented in Table 5. The exact 
values are taken from [7], 

In all these cases, it has been observed that the 
calculated values of the emergent angiilar distrilxitlon and 
the leakage are very reliable-. The calculated values 
converges to the exact values as the number of transients 
terms is increased. Thus the use of M^^s and 3^s as 
the zeros of the orthogonal polynanials giv^ by Eq* (2.20) 
and e frcm the Eg>. (2v7a) is justified. 
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Further the results will Improve if the orthogonality 
integrals (1*7) w.r*t*y. are used in obtaining the coeffi- 
cients a^^ and and as the zeros of the 

orthogonal polynomials given by Eq, (2.20) and £ from 
the Eq,. (2 •7a)* 
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Tcibl 65 2.r# 1 

Values of and'Z^ for different c's 

Q V 


G.l 

0.990542 

1.0 

8*539 

0*2 

0.981642 

1.000091 

3.9255 

0.3 

0.973524 

1.002593 

2.497 

0.4 

0.966199 

1.014586 

1.B263 

0.5 

0.959548 

1.044382 

1.4414 

0.6 

0.953451 

1.102132 

1.1925 

0.7 

0.947812 

1.206804 

1.0181 

0.8 

0,942561 

1.407634 

0.88913 

0.9 

0,937643 

1.903205 

0.7896 
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Table 2.3 




c 

e 

LEAKAGE 

exact 

PER 


OV4289 

-0.8326 

-0.8266 

-0.721 

0.2 

.0.1794 

-0.8293 

-0.8266 

-0.326 


0-.4289 

-0.7499 

-0.7427 

-0,969 

0.3 

0.1794 

» 74-44 

-0.7427 

-^.234 


0.4289 

-0.6724 

-0.6627 

-1.457 

0»4 

0.1794' 

-0.6641 

^.6627 

-0.217 


0.4289 

-0.6016 

-0.*;8B1 

-2.289 

0.5 

0^1794 

-0.5897 

-0.5881 

-0.-272 


0.4289 

-0.5363 

-0.5170 

-3.727 

0.6 

0.1794 

-0.5194 

-0.5170 

-0.46 


0.4289 

-0.4733 

-0.4466 

-5.969 

0.7 

0.1794 

-0.4495 

—0.4466 

-0.643 


0.4289 

-0,4046 

-0.3713 

-9.456 

CD 

# . 

o 

0.1794 

-0.3741 

-0.3713 

-0.753 

0.9 

0.4289 

-0.3199 

-0-2772 

-15.39 

O'. 1794 

-0.2799 

-0»2772 

. -0.957 


I 
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Table 2:.4 

^ niim pcoUaa uttfet arf 


c 

€ 

2-» J 

EXACT 

PER 

0.1 

0.4289 

-0,0079 

-0.0217 

63 . 26 ■ 

0.1794 

-0.0056 

-0,0217 

73.93 


0.4289 

-0.0539 

-0.0463 

-16,47 

0.2 

0.1794 

-0.0487 

-0,0463 

- 6*295 


0.4289 

-0.096 

-0.0745 

-29.27 

0.3 

0.1794 

-0.089 

-0.0745 

-20. 6l' 


0.4289 

-0.1362 

-0,1073 

-26.95 

0.4 

0,1794 

-0.1288 

-0.1073 

-20.08 


0-.4289 

-0,1760 

-0.1465 

-20.162 

0.5 

0-.1794 

-0.1688 

-0-.1465 

-15..^213 


0.4289 

-0.21B7 

-0.1947 

-12.31 

0.6 

0.1794 

-0.213 6 

-0.1947 

-9.699 


0.4289 

-0.2686 

-0-.2566 

- 4.691 

0.7 

0.1794 

-0,2693 

-0.2'^66 

- 4.941 


0,4289 

-0.3357 

-0.3419 

1.819 

00 

o 

0.1794 

-0.3471 

-0.3419 

-1.513 


0.4289 

-0.44919 

-0,4780 

+6.025 

0.9 

0.1794' 

-0,4758 

-0,4780 

Hh 0*4 SS 
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Table 2.5 


titf n$Xm ultli ami 



e 

2* J 

EXACT 

EER 

0,1 

0.4289 

-1.1023 

-1.087 

-1.403 

0.1794 

-1.1004 

-1.087 

-1.23 


0.4289 

-1.1827 

— Ial92 

40 ^ 7 84 

0,2 

0.1794 

-1.1891 

-1.192 

+0*242 


0,4289 

-1.2911 

-1.322 

+2.339 

0.3 

0.1794 

-1.3004 

-1.322 

+la634 


0-.4289 

-1.4396 

-1'.4880 

3.252 

0.4 

0r.l794 

-1.4519 

-1:.4880 

2.424 


0.4289 

-1.6479 

-1.7070 

3 ^465 

0,5 

0-,1794 

-1,662A 

-1.7070 

2.611 


0*4289 

-1.9533 

-2.013 

2.964 

0,6 

0,1794 

-1.966 

-2r.Ol3 

2.333 


0.4289 

-2.438 

-2, <#'^8 

1.619 

0.7 

0.1794 

-2,436 

-2.478 

1.706 


0.4289 

-3 .321 

— 3 ’• 291 

-0.930 

CD 

O 

0,1794 

-3.2646 

-3.291 

0.801 


0.4289 

-5.508 

/ 

-5.220 

-5.518 

0.9 

0.1794 

-5.242 

-5.220 

-0,427 




5is;«=.io 

SS:W3.4!) 

^s;ws,50 

4S'Ma.6Q 

N[S:i«s.-70 

??E:«a,6D 

sie:w3,9o 

'IS:«a.,40 

»S:^3,2!? 

NF5W=,30 

4S.4s.,4a, 

ME;<?a.5a 

•3S;^^s,63 

4Wa.,7D 

MS:Wa,'20 

^s:rf«s.,'30 

*IS,«I*»40 

4s:<»s,50 

Nf£:w=.,60 

4S.4S..70 

NIEWa^,€0 

MSWa-.-OO 


NfBiif*'.! 0 
!IB»* ^ 
NfEMa... 

*IE»3.,4 

»P».5_ 
SBMs,60 
»EW»,f 0 
wewa.da 


Npa.,10 

Me:*««,ao 

ra»^3o 

|i|Si(fia.,4D 


,10 

VhLs 0 
VALx 0 
VAL=-0 
\r4La-0 
VAL=-0 
VAhxmO 
VMiS-O 
VftLa*0 
V4t,a-0 
:a .20 
7 A La 0 
^TALa 0 
\rALa-0 
VALa-O 
VALa-O 
VRLa-O 
VALa-0 
VALa-0 
VALa-0 
C5s ,30 
VALa 0 
VALa 0 
VALa-0 
VALa-0 
VAL=-0 
VALa-0 
VALa-0 
VALa-0 
VALa^O 

VALs 0 
VALa 0 
VAL®*0 
VALa-0 
VALa**0 
VALa-0 
VAL»-0 
VAL=-0 
VALa-O 
C!x .50 
VALa 0 
VALa 0 
VALa-0 
VALa-0 


,64571596E-06 

.7??16532E-07 

,835379905-06 

,183977961-05 

.28434109E-05 

,380810955-05 

.471821005-05 

.556890635-05 

.636054705-05 

.535385035-05 

,707294005-06 

.662996395-05 

,146901495-04 

,227234835-04 

,304229005-04 

,376691835-04 

.444287995-04 

,507082655-04 

,184244065-04 

.232898325-05 

.227235525-04 

.500328885-04 

,771021315-04 

.102935975-03 

,127165715-03 

.149703355-03 

.170589115-03 

,442137415-04 

.447174875-05 

.561205071-04 

,121524181-03 

.185919065-03 

,247062736-03 

,304176675-03 

.357123635-03 

,406050685-03 

.879589445-04 

.570132905-05 

,116868565-03 

.247787975-03 



M£:4s,50 

'^£.^=■•60 

SS:k?=!,70 


!iSW=,10 

SPiMS 

Sl:3&*lo 

Ms:<^=.,5a 

Me.4a,'70 

^S.rfa.eO 

^£:^a.-90 


«£Ws,-10 

!^s:^^s.•20 

f«E:j^a.,'30 

'Is;^^s..40 

.'iE:Ws,60 

SIE:<#s,-70 

!«e:ws,«o 

«s;»f=,-90 

SE«».20 
»E:?#s,>3 0 
»E:»«e,40 
>JE;«a.,50 

NS:»=,70- 

'JEWa.BO 

!<E:4a..lO 
SS:«la,20 
lfw=.'30 
irEi3.,4D 
MSMS-.60 
!IE:»a,60 
S Els. 70 
S.jf 0 
MS;«a.,'90 


VAf-=-0,37S78498E-03 
Vatis-O, 49668 048E-03 
V4L=-O,6O913697E»03 
V4r.a-0.71303l5ll-&3 
„V4L=-0.80876152i-33 

V - • 0 P 

VAL.* 0,35996449E-05 
7AL*-0:22074958|-53 
VALS-0.45791240E-D3 
7ALa-O,68819765E-03 
VADa-0, 904600241-53 
7AL=-0,11050897E-02 
VALa-5,l2897078E-02 
VAtja-0. 1459351 ll-02 
7 P 

“VAfla 0.26827675E-03 
VALS-0.59061941E-05 
VALs-0,39660064E-03 
7ALS-0.80555772E-03 
VALS-0.12000921E-02 
VA^a-0.i5690713E-02 
7AL=-0,19096366E-a2 
VA£#s-0.22222855E-02 
^VAi:,s-0,25088439E-02 

gf| 

'"VAt.® 0.45594367E-03 
7Alia-0,3ll47512E-04 
7ALa-0.71ll8113E-03 
7ALa-0,l4164509E-02 
yAX<a-0.20927498E-02 
yAI^=-0,27224388E-32 
yAL=-0.33016243E-02 
VAl.a-0,38318410B-02 
yAL=-0.43166766E-02 
Ca ,90 

VkU^ 0.85529485E-03 
yAlia-0,98651250E-94 
tfALa-0.l4053942E-02 
yAliS-0,27489950E-02 
yALa-0.40301967E*02 
yALa-0,52182089E-02 
yAt.a-0,63O74524E-02 
VAli=-0,73020332E-g2 
yAl,a-0,82095449s-02 
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s program finds thie serps of 3RTS3S3MMI pOLINomi&LS 
dftSrefi N*l to 5 

IIPGICIT DOUBLE PR|?PII)« CA-H,3-2) *T,.rn.«Afv% 

OIHEflSIOM P(10O),R2Ci0O4SgKAMC153»i00),AIVCO!l0O) 

DIMENSIO?! CCl2),Ati0(12),ZEtl00)fJfl2(t23 
RE«VDC21.*)N 

FDrSIt(/I32x! 'SOLUTION OF ORTHOSOlftL POLY EQS',/) 
WRITEC24fl8'* •’ 

F3RM!1|(33X, 

REA0(2UJ)(C(iL|=|#|0) 

R£A0C21,4){RNUCI:)#I«| rlO) 

RSADai,*3CWi2Cn,I«lrlO) 

SRlIlcIlh?!) CCPilNUCI) 

F3RM4h/{25X.'Cf%,F5,2,5x.'NEW3a.%3) 
CS*C(I)/l2.D54«12Up(l.O0-0.(I))) 
5RBCS*AMU(!)fDSORr(1.00-M(I>} 

00 20 11=0.22 

ro 2i 

ftmm=<i;6D0-K!c»Asa;(i)niKii-a))/xx 

GO TO 20 




=DLOG((A»UCI)+1.00)/ANUCI3) 


DO 30 IR”}4} 

A3(iLic)=6MAivaR+in*i)+cs»AmrR+i33 

R2Ci£i)»osoRTa«a#i)) 

R2(|pCl)=IfClrICl)/R2(ia> 

00 50 


DO 50 ICsIRfll 
IFCIR.EQ.ld GO TO 51 

SIs|a=SIG«^+R2(15»|R)|R2C |S#|C) 
R2aR,lC)*(AM(IRVlC}-SIGMA)/R2CI.RyIR) 


GO TO 50 
SIGWASO.ODO 




Blnfifil 
Eraci)*. . 

0 60 3«2.S 




IR)-SIG«A) 
1 ) 


l|pc Jl»CP^CUrJ^+^*5/R2(J, J))-CR2(3-l f J)/R2CJ-lr J-i) 3 

«M!j?*R^fU+lI§+l)/R2(lr JO 
iS22^i8iTSCNtALPrBErA;,ZZ) 


S|Pl8|TSCNrALPrBErA;,ZZ) 
lo^fiio ROOTS OF P0L)fB3*?I ALS BSf SAOf -I^irERYAL HETHOD 



fK J 


3 


5 

6 


in 


u 


SUSROyTINE RO0TS(Mf.ft^|£3|rA»ZZ) 

IMPLICIT OOUBLE PREC|SIDNCA-H, 3-Z) 

gSSIfsiga 

iPli?kg!5gs“!5?a}roriH^?«‘’SSibi'’«o5it'"LSg5 

ZRsO.ODO 
00 7 1=1 fNpi 

ISO TO 3 

ZLaZb+O.OSDO __ 

ffCE.Gt.l.ODO) GO to 44 

GO' TO i ft-rvrt 

00 6 J=iiITER 


^$??¥«xxxxxxxxaixxHxuxixxxKxa. 

poiii TiPfprPPIIIt ^ 

iFCFZHALF4FZL,LE^O,0DO) GO TO > 
ZL»ZHALF 
fzl=fzhxlf 
GO TO 6 
ZR»ZHftLF 


Jtl 

'HE RODTt't# 



',EPfiETh^5 
5 GO TO 44 


dbMtl»OE. . 

^RMEjO«FPAU0 

fifexi«g§5iJi%p?ii53Sc?2iiysi4r'-"'“> 

Ifilligg 

l\ t 5=1 X-alp ( i ) ) /B|f ft Cl ) 

MOqJ) GO TO 33 

?j]-(feiLPCJ))4PCJ-l)-BS.TA(J-l)^PCJ:-2))/BETftC0) 

IFiKf 

S r ■ --f-'ij,*- 



SDWriON or ORIflOGOHM* POUf EQS' 
OSGREE OF POWIfOMIAIiB I 


o.n newob, 0,1000000004190009000+31 
HEW*. O,66787S99768501S2SOO+OO 



0,25 

MEWa- 

«<f 1^'' 

0.3 0 
M2»s> 

2!S. 

0.40 

MEMs> 



cs® 



0.7 5 
MSM®. 

:!s, 

0,80 

MEMs. 


0.90 

MSMs 


HErfOs 0,1044382934000000000+3l« 
0,6?473§589049l439540+00 

NEWOa 0,1206804258000000000+01 
L6eO45Hl7266l7801&O+OO 

NEWOs 0,1407634301000000000+01 
0.68453§323020074634D+00 

neWOB 0*|90320484|000000000+01 
|.69073§592815f 783370+00 



SOUiJIIDN OF 0RTHD50NAI, POO« EQS 
DEGREE OF POOySOMIAGR 2 


:;s 0.13 , 

M £.*!*■ 


M E < is / 


: i = 0 , 33 . 
M|»a. 
MSJ # a . 

:i= 0,43. 

MSMs . 

MEiSs. 

0.53 

mIm*. 

M !:«#*■ 


:ss 0,6 3 


M £ W «. 

MEija . 


;S5 0,73 


MEM*. 

KEM *. 


C»< 0,80 


MEM*. 

MEM*. 


:pt‘- 0,9 


ME»«. 


^0*10*500000341 00300003+31 
O,356O9§811O926418!40+OO 
0,8453864457125746430+00 

0,8458803910867: 




+0 


p ? 

932 

33556071 


00000003+01 

;3Q+00 

70+00 


NEMQa^ 0.1002592 

0,8470974652973382050+00 

MEMOs^^O, 1044382034000000000+01 
0,3622751184212170786+00 
018478624469575152050+00 

NEMO*. 0,110^32022000000000+01 
0,36459§3722526327480+0a 
0^8487800584534852550+00 

8EWgs^^O,12O6fO42|8OO0OOOOO0+Ol. 

5,371266094 
0351440595 




301000000000+01. 
ifo | 00 + 6 g 
«$604OO 


MEMO* 04l8O32Opi3OOOO0pO0+0l' 
0,1768701682198297840+00 
611537452799104357850+00 


SDtiUTION OF ORTHOGOSJiti POCif E2S 
OSGREE OF POLYMOMIACiB 3 





OOOOOOOOtDl 

OtOO 
64780f00 


0.23, NEi«jg=,,g,igoO09O8865OO3O3(JO3t-3l 
11'!* g.2n737532982853O93D+!>0 

MS|ls 0.5920490505912312090f00 
MS>J= 0,9118300872283725790+00 

:is 0.33 «E»g=^ 0,1002592888000300000+31 

«|«= 0.214614753197023037D+00 

MEitf=. 0,5929787133882230600+00 
MEWa 0.9120835909419838600+00 

0,40 ^|*^gT«,§*igi|5858l60000goOOO+3l 

MSjfis 0.2156499775384872920+00 
MEi«s 0.5940606576248683280+00 
MEMs 0,9123768413122888910+00 

0,50 NETOa, 0.1044382034000000000+31 
MEMs. 0.2168793346612831550+00 
MSMa 0,5953345562546019220+00 
«£»!=. 0.9127211368711869020+00 


G5S 0.60 NEMOS - 
MSMs 0,21835 


«BMs. 

HEWS. 


0,1102132022000000000+31 
§7102131898290+00 


S*5§6i649079Sl5^i4520+00 
0.9131350403924443550+00 


0.70 NEMOS 0,1206804258000000000+01 
HEW*. 0.22019252280333l03|D+00 
HEM*. 0,|9876?80127997|3570+00 
HEWS 0. §136520706259034360+00 


CIS 0,80 MEMOS 0,1407634301000000000+31 
HEMS 0.22258|470|588357570+g0 

M EM Si 0 , 

HEMS 0.1 


i0|2825605965l 90350+03 
)143418350i66030020+00 


0,90 NEMOS 0,1903204843000000000+01 
MEWS 0.22608§9410870984430+00 
MSm= 0.6050603175772266710+00 
MSWs, 3,9153961748201254550+00 



SDUUnON OF ORrHOGONilL POL? E3S 
DEGREE OF POOyMOMIAGBs 4 


34401580618392472390 + 00 
0.4170532671527 


0.7235 

0.9429 


110629! 

1912831 


^325880+00 
37I3208D+00 
86555120+00 


::s 0,10 
MSjrfs 
MSWa 
MS.Wa 
MS#ls. 

CIS. 0.20 

MSMs 
MEa’s 
MEWs. 

MBW* 

C^s 0.3 0 

MBM=- 
MSitfs 
MSMa. 

MSW=. 

C^ 0,40 
Hsrrfs 

ME«5 
MSWs 
MEMs, 

C^ 0.60 

MSMa. 

MEM*- 
«SW*. 
mux. 

CIS 0.60 

Pm 

«flfs 0l72?Oi0567663O6S2&5D+00 
MSMs 0.9438041710283869210+00 

Cl= 0.70 MEMO® 0.1206804258000000000+01 
J4SMS 0,1446358850942488060+00 
^EMs 0,4240065454014256840+00 
HSWs 0.7281708432401501360+00 
HSMs 0,9440778142943600000+00 


NEi^Os^ 0.1000090886530000000 + 01 
0.1406199551962345140+00 
0.4177901417566317840+00 
0,7240835772219725200+00 
0,9431158498136937880+00 

NEMOS0.1002592888000000003+31- 

0.1411633422289014560+00 

0.4186447407100786230+00 

0.7246494792434532450+00 

0.9432493593169056140+00 

»E«Os 0,1014585816000000000+01 
0.1418069446390290980+00 
0.4196446914644184300+00 
0.7253067616959015160+00 
0,9434039015959569980+00 

KSMO* 0.1044382034000000000+01 
0.14257§88576081 91900+00 
0.4208257379290444080+00 
0,726Q8O263327O9414§D+g0 
0.9435855242754769310+00 



CJs 0.80 HEi^Os 0.1407634301000000000+01 

MSMs 0.1461 iHi 33027038310+00 
MSMs 0.4263207460109697430+00 
MSWs 0,7297124860084295520+00 
MBMs 0.9444438452246686240+00 

:= Q.90 ?«SWO= 0,1903204843000000000+01 

MEMs i),14826§279i88025394D+00 
MSMs 0.4297582142960891370+00 
MSM= 0,7320499983074114440+00 
MEMa 0,9450056097553897450+00 



SDbUrlOM OF ORTHOGONAL PObT EOS 
DEGREE OF POLTSOmALfe 5 


:iB o.iD 

MSWs 

MS.»s- 

ME)#*- 

MEWS. 

«E^=- 

0,2 3 
MSMs. 
MS.vfa- 
ME-Ns 
MEW*- 
MEM=. 

0,30. 

MSMs- 

MSiia. 

MENS. 

MENS, 

MENS. 


^}m 

1^' 


0.40 

MENS, 

MENS. 

MSls^ 

MEMs. 

ME.MS. 


C)s. 0.50 , 
MENS. 
MEW*. 
MEW*. 
MBW». 
MBW* 

:?s 0.60 

MEWS. 

MEWS. 

MEW* 

MSW= 

MEWS 

:Jc. 0,7 0 

MEW*. 

MEWS. 

MEW*. 

MEWS 

MEW* 

>^hm 'A A 

MEWS. 

MEWs- 

MEWS. 
MEW* 
' MEWS. 


0.9879055168295407220-01 
0,3050450900784198900^00 
0.5625175393714016560+00 
0,8022611808893998420+00 
0.960250822436864841 0+00 

MEWOs 0,1000090836500000000+01 
0.99086l4025895680510-]0r 
0.3056308083756448480+00 
0.5630786717007140410+00 
0.8025723662616655960+00 
0,9603193882565392450+00 

NEWOs 0.1002592888000000000+01 
a,99435§4445340284790-01 
0.3063133762367215250+00 
9,5637264486536878410+00 
0.8029294870797457430+00 
0.9603978425893728890+00 

NEWOs 0. 1014585816000000000+01 
0.99850^07 7345881460-01 
0.3071152116794110050+00 
0.5644816674070170870+00 
0,8033439721035392760+00 
0.9604887006407807390+00 

NEWOs 0,1044382034000000000+01 
9.|0034l756833874271D+00 
0.3080643573186623760+00 
0,5653721746904921040+00 


0.8038317968657793280+00 

0.9605955539616843500+00 


0 ^ 11 ^ 132022000000000+01 


:?= 0 , 


,90 

MS.Ws. 

ME-Ws. 

MSMs. 

MEWS. 

MSWs. 


NEMO* . 

0.1009397390964295480+00 
0,3092053973570727980+00 
0,5664431367436918680+00 
0,8044191307937580890+00 
0,9607243112066498770+00 

NEWOs 0.1206804258000000000+01 
0,1016738744041504110+00 
0,3106180138132913280+00 
0.5677757877881958870+00 
0.8051532139670598580+00 
0.9608856593582459070+00 

NEWOs 0,1407634301000000000+01 
0,1026244643628160700+00 
0,3124649167704774310+00 
0.5695374744656874110+00 
0,8061320436750975230+00 
0.9611018603263801230+00 

NEWOs 0,1903204843000000000+01 
0.10399^0121848707790+00 
0,|l51821243493pifSD+00 
0.5721829549849644550+00 
0,8076250707217695920+00 
0.9614345618316292570+00 
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